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Abstract. In this article we prove solvability results for L 2 boundary value prob- 
lems of some elliptic systems Lu = on the upper half-space R™ +1 ,n > 1, with 
transversally independent coefficients. We use the first order formalism introduced 
£SJ ' by Auscher-Axelsson-Mclntosh and further developed with a better understanding 

of the classes of solutions in the subsequent work of Auscher-Axclsson. The inter- 
esting fact is that we prove only half of the Rellich boundary inequality without 
knowing the other half. 

O 

Q 

t^- 1. Introduction 

The goal of this paper is to study Rellich type estimates for some elliptic systems 
^ ■ by using the first order formalism introduced in [3J. 

. We begin by recalling the classical Lax-Milgram existence theorems of variational 

(or energy) solutions for divergence form second order partial differential equations 
and make clear how we understand them in unbounded domains. It is well-known 
that both Neumann and regularity problems have unique solutions in the energy 
class. We reformulate this using the DB formalism in order to represent these solu- 
tions via a semi-group. We then obtain factorisations of the Dirichlet to Neumann 
map and of its inverse in H~ 1 / 2 topologies which are related to (abstract) bound- 
ary layer potentials. This part is completely general and true for any such elliptic 
^ . system. 

^C) \ In the last part, we assume that the coefficients are block triangular (see below 

; for definition). Although such an assumption is rather restrictive, it brings up a 

phenomenon which we think interesting on its own. Using these factorisations, we 
prove that, in the triangular situation which corresponds to making the conormal 
derivative proportional to the transversal derivative, the Neumann to Dirichlet map 
is bounded in the L 2 topology, that is, solutions (we shall make clear the meaning 
of solutions) satisfy half of the boundary Rellich estimate 

||V tan n|| 2 < C\\d A u\\ 2 

which implies that the Neumann problem with L 2 data is solvable. For the adjoint 
situation, it is the Dirichlet to Neumann map that is bounded on L 2 , hence the 
opposite Rellich inequality for solutions 

IIMh < c\\v tanU \\ 2 

holds so that the regularity problem with L 2 data is solvable. In each case, we do 
not know about boundedness of the inverse, that is, we do not know the other half 
of the Rellich estimate (and we do not expect it by any means) which is the alluded 
to phenomenon: usually Rellich estimates are proved both ways. We also show that 
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the Dirichlet problem with L 2 data is solvable in the same block triangular situation 
as for the Neumann problem. 

We remark that the conjunction of the two situations is the block-diagonal case 
for which 

||V tan n|| 2 » \\d A u\\ 2 

is known to be equivalent to the Kato square root estimate [I] . 

As the reader might notice, a feature of our proof is the use of one equivalent 
formulation of the Kato square root estimate (for an auxiliary operator) which shows 
a tight connection between Rellich estimates and Tb technology. 

In a subsequent paper of two of us (the first and third named author) with S. 
Hofmann, Neumann solvability results for Hardy space data are obtained in the 
triangular situations of this paper in the case of systems of such pde's satisfying the 
De Giorgi regularity condition. 
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Mcintosh acknowledges support from the Australian Government through the Aus- 
tralian Research Council. Mourgoglou was supported by the Fondation Mathemati- 
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2. NOTATION 

The system of equations is 



n m 

(Lunt,x) = J2 ^^(^(^/(^x)) = 0, 



= }^2^ d 

i,j=0 13=1 

in M^ +n = (0, oo) x R n , where d = §- t and 9; = if % — 1, . . . , n. We assume 

(2) a = {A-f{x))tf=tT e L °o( K "; £( c2m ))' 

and that A is strictly accretive on the subspace "H of L 2 (R n ; C( 1+n ) m ) defined by 
{fj l )i=x,..., n is curl free in W 1 for all a, that is there exists A > such that for all 

fen ' 

(3) 



The system ([T]) is always considered in the sense of distributions with weak solu- 
tions, that is, H\ oc solutions in the particular domain under consideration. 

We stress that our results are valid for any m but set for notational simplicity 
m = 1 from now on. In this case, the accretivity condition above is equivalent to 
the usual pointwise accretivity condition 

n n 
i,j=0 i=0 

It is convenient to write A in the block form 

a(x) b(x) 



A(x) 



c(x) d(x) 
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where a is scalar- valued, b, c vector- valued and d n x n matrix-valued. Call A the 
set of 2 x 2 block matrices A with these properties. 

All our estimates in this article depend only on ellipticity constants A = ||^4||oo 
and the largest A in the accretivity inequality (EJ). For non negative quantities a, b, 
the notation a < 6 means a < Cb where C is a constant depending on the parameters 
at hand. The notation a m b means both a < b and b < a. 

3. Energy solutions 

In this section we recall the usual construction of variational or energy solutions. 
Although this is fairly classical and mostly follows [3J, we need to stress a few points 
and set up some notation. This uses the homogeneous Sobolev space if 1 (M^ +n ) of 
u G L / 2 oc (lR+ +n ) such that V t , x u G L 2 (IR+ +n ), equipped with the semi-norm IHI 2 ^ : = 
JJ R i +n | Vt, x u\ 2 dtdx (note that we could have supposed u G X>'(IR+ +n ) as it is not 

hard to check that u can be identified with a L 2 oc (IR^ +n ) function when Vu G L 2 ). 
This is a Banach space modulo constants. 



Lemma 3.1. Co°(M+ +n ) is dense in if 1 (M+ +n ). The trace on R n is bounded from 
H 1 (R ] ^ rn ) onto the homogoneous Sobolev space if 1 / 2 (IR n ) ; which we define as the 

closure ofC yo (M n ) for the semi-norm [J Rn (with f designating the 

Fourier transform) . 



Proof. Once the density is shown, the trace result is classical. Let u G i/ 1 ) 
As for any compact subset K of W 1 and < to < ti < oo, 

\u(t\, x) — u(to, x)\ 2 dx < (ti — to) // \d t u(t, x)\ 2 dtdx, 

. ./ICP 1 + n 



A 



u extends to an element u G L 2 oc (M^_ +n ). By the reflection principle, we can extend 
u to an element in if 1 (M 1+n ) which we still call u. We claim we can find a sequence 
Uk of L 2 (E} +n ) functions with compact support that converges to u in ij 1 (R 1 f +n ). 
Indeed, by Mazur's lemma, it suffices to have a sequence with weak convergence, 
that is such that Vn^, — > Vu weakly in L 2 . It is easy to see using Poincare inequality, 
and this is where we need to know a priori that u G L 2 oc , that = (u — c^)^ has 
this property when is the mean of u on the ball -8(0, 2 fc+1 ) and (fk{x) = ip(2~ k x) 
with (f is a smooth function that vanishes outside -8(0, 2) and that is 1 on -8(0, 1). 
Finally, it suffices to mollify this sequence to conclude. □ 

Theorem 3.2. Given £ G ij _1 / 2 (M n ) ; there exists u G .8 1 (]R^ +n ) ; unique up to a 
constant, such that ff R i+n AVu ■ V(pdtdx = (£, ip), for any <p G -8 1 (R^ t +ri ) with trace 

(p. This function u is a weak solution of Lu = in M^_ +n . We define the conormal 
derivative of u at the boundary to be d VA u\t=o = —£ (we use the inward unit normal 
convention) . 

We say that u is the energy solution of the Neumann problem for Lu = with 
Neumann data —i. 

Proof. This is just the Lax-Milgram theorem in the Hilbert space H 1 (IRj l f +n )/C using 
that H^/ 2 (R n ) is the dual space of H^ 2 (R n ). □ 
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Theorem 3.3. Given f E H 1 ^ 2 (W a ), there exists v E H 1 (M}^~ n ), unique up to 
a constant, such that Lv = in M^_ +n and V\ t= o = f with equality in ij 1 / 2 (IR n ). 
Furthermore, there exists £ E ij _1 / 2 (R n ) such that JJ m i+n AVu-V(p dtdx = (£, cp) for 

any p E i/ 1 / 2 (IR n ) and any extension <p E H l (M}+ n ) of p. 

We say that v is the energy solution for the regularity problem Lv = with data 
V x f and we have d„ A v\ t =o = —£■ 

Proof. Pick an extension w E ^(R^ 71 ) of /. Let i^(R+ +n ) be the subspace of 
if 1 (R^. +Tl ) consisting of all u with constant trace on R n (alternately this is the closure 
of Co°(]R^ +n ) in if ^R^ - ™)). By the Lax-Milgram theorem, there exists a unique 
u E ^o( K + + ") solving 



AVu ■ V<p dtdx = -II AVw ■ V<j> dtdx 

l+n J J^}+n 

for all E ijo(R+ +n )- Then v = u + w is the solution. 

Next, the integral JJ R i+« AVu ■ V(f>dtdx depends only on the trace modulo con- 
stants of E if 1 (R+ +Tl ). Thus, the map AVu ■ V0 dtdx is bounded from 
ij 1 / 2 (R rt ) to C and this defines I. □ 

Observe that V x is injective with closed range from ii" 1 / 2 (R n ) into ij" 1 / 2 (R n ; C"). 
We set "Hy the range of this map and for a reason that will become clear later also 
set iiT^ 2 = if _1 / 2 (R n ). By a Fourier transform argument, one sees that "H^ 2 = 

n(Hl 1/2 ), where TZ = V(-A)^ 1 / 2 is the array of Riesz transforms on R n (the Hilbert 
transform if n = 1) and A is the ordinary self-adjoint Laplace operator on L 2 (R n ). 
With this notation, one defines the Neumann to Dirichlet map 

(4) r A ND i = v x u lt=Q , i e n~ 1/2 

where u is the energy solution of the Neumann problem for Lu = with Neumann 
data — i and the Dirichlet to Neumann map 

(5) Tj) N g = d VA V\ t = Q , g E U^, 

where v is the energy solution of Lv = of the regularity problem with data g. 

Theorem 3.4. r^ D is a bounded and invertible map from r H ± 1 ^ 2 onto T-L^ 1 ^ 2 with 
inverse Tp N . 

The proof is an obvious consequence of the above results with our definitions. We 
have ^N D (d,y A U\t=o) = V x u\ t=0 and T^ N (V x u\ t=0 ) = d UA u\ t=0 for any energy solution 
u in the upper half-space of Lu = 0. 

We finish this section with a standard result. 

Lemma 3.5. Let u E F 1 (R 1+n ) be a solution of Lu = in R 1+n . Then u = 
(modulo constants). 

Proof. By definition, we have 

AVu-V~4>dtdx = 
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for all <p G C£°(R 1+n ), hence for all <p G H 1 ^ 1 *") by density as in Lemma EH We 
conclude taking <p = u and using the accretivity of A. □ 

4. The first order formalism 

Following [3] and [TJ , we can characterize weak solutions u to the divergence form 
equation (JTJ), by replacing u by its conormal gradient Vam as the unknown function. 
More precisely (CD) for u is replaced by ((Zj) for 

d VA u(t,x) 



F(t,x) = V A u(t,x) 



V x u(t,x) 



where d VA u 



(AVt,x u )±: that is the first component of AV t x u. This is the inward 
conormal derivative of u for the upper half-space and the outward conormal deriva- 



tor vectors in C 1+r \ 



tive for the lower half-space. Here we use the notation v = 

where v ± G C is called the scalar part and v n G C n the tangential part of v. For 
example, d t u = (Vt, x u )j- an d V x m = (Vt,x u )\\- 

Proposition 4.1. The pointwise transformation 



(6) 



A\-> A:-- 



"1 


0" 




a 


b 


-1 


a 1 —a x b 


c 


d 







1 




ca -1 d — ca~ 1 b 



is a self-inverse bijective transformation of the set of matrices in A. 

For a pair of coefficient matrices A = B and B = A, the pointwise map V^it i— > 

(RF) 

F = V au gives a one-one correspondence, with inverse F i— y Vt x u = r 

L 

tween gradients of weak solutions u G if 1 Q C (Mi 1 ( +ri ) to (QJj and solutions F G Lf oc ( 



be- 



jl+n. (nl+n\ 



of the generalized Cauchy-Riemann equations 

div, 



(7) 



d t F + 







where the derivatives are taken in the 



BF = 0, curL.Fi, = 0, 
i +n distributional sense. 



This transformation was introduced in [3] and the proposition is proved in this 
generality in pQ. We shortly review the L 2 theory in [1]. Denote by D the self-adjoint 
operator on H = L 2 (W l ; C 1+n ) defined by 



D := 





-v. 



diy, 




D(D] 



D(V) 
D(div) 



The closure of the range of D is the set of F G % such that curl^i = 0, that 
is R(D) = "H . It is shown in [7] that the operators DB and BY) with respective 
domains 5 _1 D(D) and D(D) and ranges R(D) and BR(D) are bi-sectorial operators 
with bounded holomorphic functional calculi on the closure of their range "H and 
B^ lr H°. Observe the similarity relation 

(8) B(DB) = (B~D)B on D(DB) 

that allows to transfer functional properties between DB and BY). In particular, if 
sgn(z) = 1 for Re 2; > and —1 for Rez < 0, the operators sgn(DI?) and sgn(i?D) 
are well-defined bounded involutions on "H and B^T-L respectively. One defines 
the spectral spaces "H^b = N(sgn(D£>) =1= I) and = N(sgn(_BD) =1= I). They 

topologically split "H and B~ lr H° respectively. The restriction of Y)B to the invariant 
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space is sectorial of type less than 7r/2, hence it generates an analytic semi- 

group e~ tr>B , t > 0, on it. Similarly, the restriction of BY) to the invariant space 
"H^d is sectorial of type less than n/2, hence it generates an analytic semi-group 
e" tBD , t > 0, on H i&. 

Theorem 4.2. Let u G H} oc (E}+ n ). Then, 

(1) u is a weak solution of Lu = with ||iV(Vu)||2 < oo if and only if there 
exists Fq G "Hdb such that Va« = e~ tDB F . Moreover, F is unique and 
\\F Q \\ 2 n\\N(Vu)\\ 2 . 

(2) u is a weak solution of Lu = with JJ R i+n t\Vt, x u\ 2 dtdx < oo if and only if 

there exists F G such that V au = De~ tBr> F . Moreover, F is unique, 
\\F0W2 ~ {jj^+nt\V t^u^dtdx) 1 / 2 and u is given by u = — (e~ tBT> F ) ± + c for 
some constant c G C. 

The if part was obtained in [3] and the only if part in [TJ Theorems 8.2 and 9.2]. 
Here N(g) is the Kenig-Pipher modified non-tangential function || A^Vw) H2, where 

N(g)(x) := supt-( 1+n » 2 \\f\\ L2mt>x)) , x € W 1 , 
t>o 

with W(t,x) := (c^t, c t) x B(x;cit), for some fixed constants c > 1, C\ > 0. 

Remark 4.3. Although we do not need that, the same proof shows when coefficients 
are i-independent that for the equivalence of (1) to hold one could replace ||iV(Vu)||2 
by the weaker condition sup <>0 (| J t 2 * || V SjX m||2 ds) 1 ^ 2 or by sup t>0 ||Vt )a; u||2 or even by 
the square function (/Li+n t\d t Vt,x u \ 2 dtdx) l l 2 , so that in the end all these quantities 
are a priori equivalent for weak solutions. 

In (2), the function Fq is formally built as DFo = Vyi^|t=o (which belongs to an 
adapted Sobolev space of order -1: we shall make this more precise). 

The spectral spaces with negative signs correspond to estimates for solutions to 
Lu = in the lower-half space, and the similar statement holds using the semigroups 
e mB and e tBD . 

Our aim is to extend this formalism to Sobolev spaces. However, a difficulty 
is that R(BD) is really an L 2 object as it depends on B. We shall modify the 
setup to prepare this extension. Recall that "H = R(D) is a closed subspace of 
U = L 2 (R n ; C 1+n ). Let S = D| W o with domain D(D) n H°. Then S is a one-one, 
self-adjoint operator. Define n as the orthogonal projection from "H onto "H . It 
is the identity if n = 1 but not otherwise. Let B be the operator on "H defined 
by Bu = UBu = UBUu,u G "H . As B is a strictly accretive operator on "H (for 
equations this is true on % but only on T-i° for systems, that is, when m > 1), the 
restriction of n on B7i is an isomorphism onto "H and B is a strictly accretive 
operator on %°. 

Define 

T:U° ^H°, T = BS = TlB~D\ n o with D(T) = D(S) 

and 

T:U° ^U°, T = SB = Dn£| H o = DB ]H o with D(T) = B'^S). 
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Theorem 4.4. (1) T is a one-one bi-sectorial operator with bounded holomor- 
phic functional calculus on H°. A function u £ if/ oc (M^_ +n ) is a weak solution 
of Lu = with \\N(Vu)\\2 < oo if and only if there exists Hq £ 1-L^ + such 

that Vau = g^—Hq. Moreover, Hq is unique and H-H0II2 ~ || iV( Vti) || 2 - 
(2) T is a one-one bi-sectorial operator with bounded holomorphic functional 
calculus on Tt°. A function u £ Hj- 0C (Rr^" n ) is a weak solution of Lu = 

with J f R i+ n t\ \ 7 tx u\ 2 dtdx < 00 if and only if there exists Hq £ "H^' + such that 
V au = Se~ tT H . Moreover, Hq is unique, \\Hq\\ 2 ~ *l V t , x u\ 2 dtdx) l l 2 

and u = — (e~ tT H ) ± + c for some constant c. 

Proof. We begin with the first point. As T is the restriction of DB to "H , it is 
a one-one bi-sectorial operator with bounded holomorphic functional calculus. In 
particular, the spectral spaces coincide with H^s- ^ ne function H is nothing 
but Fo in Theorem 14.21 (1). This proves the first point. 

Let us turn to the second point. That T is one-one, bi-sectorial with bounded 
holomorphic functional calculus on "H follows from the similarity relation BT_ = TB. 
To obtain the new equations for u is not as direct. Denote by U : BH° — > 7i° the 
restriction of IT on BU°. For G £ H°, we have UiU^G - G) = so U~ l G £ D(D) 
if and only if G £ D(D). A calculation shows that TG = UBBG = UBDU^G. So 
T is also similar to the restriction of I?D on B7i°. Next, the equation F = e~ tBr> FQ 
with F £ is equivalent to UF = e~ tT UF with UF £ . Hence, 

Se- tT UF = T>Ue- tBD F = De~ tBr> F 

and 

(e- tT UF ) ± = (Ue- tB »F ) ± = (e"^ D F ) x 

as U leaves the scalar component invariant. This concludes the proof with Hq = 
UF = HF , where F is the function specified in Theorem |4.2| (2). □ 

The point is that both operators T and T_ act on the same space "H which is 
independent of the action of B. Thus, we use the relation between quadratic esti- 
mates and interpolation developed in Section 8] for bi-sectorial operators and we 
strongly suggest the reader to have this reference handy from now on. 

For s £ R define W as the completion of "H for the quadratic norm 

where ip is a suitable holomorphic function on bi-sectors, for example ip(z) = 
z k e~ z sgn ( 2 ) , Re z 7^ and N 3 k > max(s,0). Remark that from the spectral theo- 
rem U° = U° and it can be checked that \\F\\ S , S = c^ s \\\S\ s F\\ 2 where \S\ = (S 2 ) l/2 . 
Note that S extends to an isomorphism between 7-L s and "H 5-1 . Classically, the 
intersection of lri s for s in a bounded interval is dense in each of them. 

We define similarly for lri s T and lri s T replacing S by T and T. The quadratic 
norms are equivalent under changes of suitable i[) that are non degenerate on both 
components of the bi-sectors. 

Note that \S\ preserves the normal and tangential components so we can write 
"H s = Tit* © "Hf which agrees with our earlier notation when s = and s = —1/2. 
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If n — 1, we have S 



D 





-oL 



4 
o 



so the quadratic norm defines the usual 



homogeneous Sobolev space if s (M;C 2 ) which is also the completion of D(|5| s ) for 
the homogeneous norm |||S| S F|| 2 . It follows that H s ± = "Hjj = H S (R). 

For n > 2, recall that 1Z = V(— A) -1 / 2 denotes the array of Riesz transforms and 
let 1Z* = — (— A) _1 / 2 div be its adjoint. The operator 



V = 

is an isometry with inverse 

V~ l = V* 





-n 



: L 2 (R n ;C 2 ) ->■ H c 





-n* 



:H° ^L 2 (M n ;C 2 ) 



and 



Thus, U\ 



II. A calculation using V 
V^SV = 



K(-A) 1 / 2 shows that 





-A) 1/2 



(-A)V* 




and ?tf = ft£P 



where H s 



is the usual homogeneous 



Sobolev space with semi-norm ||(— A) s / 2 /|| 



t|T|^||2 dt j 1 / 2 



Proposition 4.5. (1) For s G M., /or a// bounded holomorphic functions b in 
appropriate bi-sectors, b(T) extends to a bounded operator on l-i s T . In partic- 
ular, this holds for sgn(T) which is a bounded self-inverse operator on lri s T , 
and T and \T\ = sgn(T)T extend to isomorphisms between l-i s T and l-Lff 1 . 
The operator \T\ extends to a sectorial operator on H^. 

(2) % s topologically splits as the sum of the two spectral closed subspaces = 
A/(sgn(T) - J) and Wf = A/(sgn(T) + J). 

(3) The same two items hold with T replaced by T. 

(4) For < s < 1, l-i s T = T-i s and for — 1 < s < 0, l-i s T = W with equivalence of 
norms. 

(5) Furthermore, for —1 < s < 0, we have for \\F\\s, s ~ {Jq 

Proof. The first four items are contained in [HI Theorem 8.3] and the previous sec- 
tions therein with the exception of the cases s = 1 and s — — 1 of (4). 

For s — 1, observe that ||-F||t,i = c ll 1^1-^11- Thus, the bounded holomorphic 
functional calculus of T on "H implies that |||T|F|| m \\TF\\ (SJ Theorem 8.2]. 
Finally \\TF\\ « \\SF\\ « \\F\\ s ,i. 

For s = — 1, we observe the intertwining relation T_S = ST (as unbounded op- 
erators) so the functional calculi also intertwine. Hence, ^ t (T)SF = Sip t (T)F = 
t~ l B- l i) t (T)F for F G D(T) = D(S) with $(z) = zijj(z). This implies that 
\\SF\\t,~i « ||-F||r,o ~ \\ F \\s,o ~ H-S^lk-i and the result follows from the fact 
that the functions SF, F G D(T) = D(S'), form a dense subset of *H . 

To prove (5), we proceed as above and a calculation with ip(z) = ze~ zsgn ^ shows 
that e-^SF = Se^F = t^B^ipt^F for F G D(T). It follows easily for 
appropriate F and using < s + 1 < 1 that 



r 2s ||e 



dt 



1/2 



IS^Fl 



T,s+1 



IS^Fl 



S,s+1 
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□ 

We can give another useful characterization of solutions with square function 
estimates in terms of the negative Sobolev space. 

Corollary 4.6. Let u G Hl oc (M}f n ). Then u is a weak solution of Lu = with 
jj^+ntl^t^u^dtdx < oo if and only if there exists Hq G ' + such that Vau = 

e~*-iJ - Moreover, Hq is unique and ||if ||T,-i ~ {j J R i+nt\V t ,x u \ 2 dtdx) 1 / 2 . 

Proof. This is a reformulation of the previous results. By Theorem 14 .44 (2), given 
u, we have Vau = Se~ tJ Hq = e~—SHo, where e~ l — is now the semi-group extended 
to Setting Hq = SH , we have H G T-L~ T l,+ as Hq G and ||i?o||T,-i ~ 

II #0 115,-1 ~ ||-f^o II 2 by Proposition 14.51 (4). 

Conversely, let H G H^' + and set H = S^H G H ^. Then e^H = 
e^—SHo = Se~ tT H is the conormal gradient of a solution u by Theorem l4.4l (2). □ 

Now, we prove a representation for energy solutions 

Proposition 4.7. Let u G if J 1 0C (R+ + "). Then u is a weak solution of Lu = with 
jj R i+n \W t, x u\ 2 dtdx < oo (i.e., u is an energy solution) if and only if there exists 

H G 'H t 1 / 2 ' + such that Vau = e - *— Hq. Moreover, Hq is unique and ||-£fo||r,-i/2 ~ 
(JJ R i+n \Wt,xu\ 2 dtdxf/ 2 . 

Proof. Let us prove the converse first. If Hq G "H^ 1 ^ 2 ' -1- , then H £ = e~ 6 -H G 
'H T 1 ^ 2 ' + n "H^/ 2 ' + C for any e > 0. By Theorem 14.41 e~ l —H £ is the conormal 

" — 1 /2 -f~ 

gradient of some solution u £ . As H e converges to H in H, T ' , it is easy to conclude 
that u £ converges to a solution u is the energy space and (ff R i+n \ V t^u^dtdx) 1 ^ 2 ~ 

||-^o||r,-i/2- 

Assume now that u is a weak solution of Lu = with ff R i+n \Vt, x u\ 2 dtdx < oo. 

Set H = VA^|t=o ^ "H -1 / 2 - Decomposing H = Hq + Hq according to the spectral 
decomposition "H -1 / 2 = - H T 1 ^ 2,+ © 1-L T 1 ^ 2, from Proposition 14.51 (2) and (4), and 
using the implication just proved, Hq is the trace of the conormal gradient of an 
energy solution u + in the upper half-space and, by the same result in the lower 
half-space, Hq is the trace of the conormal gradient of an energy solution u~ in the 
lower half-space. It is then easy to see that the function v, defined by v = u — u + 
on the upper half-space and v = u~ in the lower half-space, is an energy solution 
of Lv = in M 1+n . By Lemma [3.51 v = (modulo constants) so that u = u + and 
Hq = Hq. Thus, Vau — VyiM + = e~*— Hq from the first part of the proof. □ 

/ 



Corollary 4.8. The elements of H. T ' are 



be written 



T 



, f G H^^ 2 . They can also 



g ' 9 t ri\\ 

Proof. It suffices to check the first representation by Theorem 13.41 By the previous 

* — 1/2 4- 

result, Hq G H, t ' if and only if Hq = W AU\t=o for some solution u in the energy 

' — 1/2 

space. So if / is the conormal derivative dAU\ t =Q in 'H ± then the tangential gradient 
at t = is T^ D f by ©. 
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— 1/2 

Conversely, if / G T-t± , then we let u be the energy solution to Lu = with 
Neumann datum /. We know Tj^ D f = \7 x u\ t =o and that V^w = e - *— H for some 

f 



• 1/2-4- 

Ho G T-L T ' . It follows that 



H G 'H T 1/2,+ . □ 



Remark 4.9. It is interesting to compare Theorem 14.44 Corollary 14.61 and Propo- 
sition 14.71 The use of T allows the same representation for conormal gradients of 
solutions and the only difference is the space to which the trace belongs. Solving 
Neumann, Dirichlet and/or regularity problems amount to proving boundedness of 
Neumann to Dirichlet map and/or the Dirichlet to Neumann map with different 
topologies on the boundary. 

Remark 4.10. In [3J section 5], some a posteriori identification of the solutions 
coming from the DB formalism in L 2 is made with energy solutions, upon some 
further assumptions. But note that energy solutions require a different trace space. 
Here, the extension of this formalism to Sobolev spaces allows to prove a priori 
representation (Proposition 14.7ft for energy solutions. 

Remark 4.11. We note that Proposition 14.71 was stated without proof for systems 
on the ball in [2], even with some radially dependent coefficients. While this article 
was in its final stage, we learned that this kind of representations of solutions with 
data in Sobolev spaces was pursued in more generality by Andreas Rosen [12], and 
has some overlap with ours. 

5. The operator theoretic lemma 

The main ingredient in our proof is the next lemma which will provide a factorisa- 
tion of the boundary maps with simpler operators to analyze. The following lemma 
is essentially taken from [HI Section 6]. 

Lemma 5.1. Let Xi, X2 be Banach spaces and let Z = X\ ©X2 whose elements are 

G C(Z) such that 



written as 


u ± 


. LetS = 


$11 $12 




My 




$21 $22 


(1) S 2 - 


= Iz 







(2) There exists c G (0, 1) such that for all u G N(S ± I), one has 

C^ 1 ||M x || Xl < ||W||||x 2 < c||ul|Ui- 

Then, S12, S21, s n S22±-f o^e one-one with closed range in the respective C(Xi, Xj). 

Proof. Let P ± = |(/± S) be the pair of projectors on Z associated to S by (1) and 
Q± the pair of projectors on X\ © {0} and {0} © X 2 . The assumption (2) implies 

||Q+P ± u|| w ||Q-P ± u|| w HP^H, u G Z. 

We infer that 

\\P + Q±u\\ « \\P-Q ± u\\ « ||Q ± w||. 

For example, using Q-Q+ = 0, then < ||P + (5+w|| + ||P"Q+m|| < ||P + (5+m||H 

\\Q_P-Q + u\\ < \\P + Q + u\\ + \\Q^P + Q + u\\ < \\P + Q + u\\ < \\Q + u\\. 
Going further, one easily sees that 

\\Q-P + Q±u\\ w \\Q+P + Q ± u\\\ w \\Q±u\\ w \\Q_P-Q±u\\ ps \\Q+P~Q ± u\\. 
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Thus, let G X 2 and u 







Q^u. One checks that 



Q+Su = Q+(2P + - I)u = 2Q + P + Q„u 



and 



hence 





(/ ± s 22 ] 



2Q_P ± u = 2Q^P + Q_u, 



\s 12( f>\\ Xl tt\\(I±S 22 )(f>\\ 



x 2 - 



Similarly, with G X\ and u 







<5+m, one obtains 



|s210|k«||(/±SllMk 



□ 



Lemma 5.2. 



'1) T7ie operator sgn(T) satisfies the hypotheses of the above lemma 



on l-i T 1 ^ 2 = % 1 / 2 = i-i^ 2 ©?^|| ^ 2 . Moreover, the operators S\ 2 {T), s 2 i(T), 
s ii{H) =t ^; S 22(X) =t I are invertible. 
(2) TTie same conclusion holds replacing T by T on "H^/ 2 = "H 1 / 2 = "H^ 2 © "HZ 2 - 

Proof. We prove (1). Recall first that U^ 12 = n~ 1/2 = HI 1 ' 2 © H~ 1/2 follows from 
Proposition 14.51 (2) and (4). Next, equality sgn(T) 2 = / on T-i^ 1 ^ 2 holds by the 
bounded holomorphic functional calculus and sgn 2 (z) = 1 on C \ R. Recall that 



N(sgn(T) - /) = Ut 1/2,+ and by Corollary EH] and setting T+(T) = Tf fD : U 
"Hfi" 1 ^ 2 ! , one has for any u G "H^ 2 



7-1/2 



— >■ 



G N(sgn(T) -/) = ?£ 



-1/2, H 



Similarly, the Neumann to Dirichlet map for the lower half-space yields the corre- 
sponding operator T~(T) characterized by for any u G "H^ 2 



tin = r (r_)u A 



G N(sgn(T) + J) = ft 



-1/2, 



Details are left to the reader. Thus, the second assumption of Lemma \5. II is granted 
and we have lower bounds for all six operators in the statement. Thus, when A = Id, 

" H 



T = S and so sgn(T) = sgn(S) 



and sgn(T) = sgn(S') 



H* 

-n 

-K 



if n = 1 where H is the Hilbert transform 



on "H if n > 2. Thus, the six operators are 



invertible in this case and in particular onto. The ontoness for general A follows 
from the method of continuity. This proves (1). 

To prove (2), we observe that the intertwining relation 5*sgn(T) = sgn(T_)S ex- 
tends to "H^/ 2 = "H 1 / 2 using that S is an isomorphism from "H 1 / 2 onto "H" 1 / 2 . Thus, 
the conclusion follows straightforwardly from (1). □ 

Lemma 5.3. In C(H~ 1/2 , H~ 1/2 ), 

T A ND = s 12 (T)~\l - s u (T)) = (I - s 22 (T))- l s 21 (T). 
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Proof. We have seen that the operators r ± (T) are bounded and invertible. To obtain 
the formulas, one uses their operator defining relations 



's n (T) s 12 (T) 
S2i(T) s 22 (T) 



± 



7 

r ± (r) 



i 

r ± (T) 

and then solve for r ± (T) using the invert ibilities of the operators in Lemma I5T21 We 
conclude using Tff D = T + (T). (As mentioned, the operator T~(T) is the Neumann to 
Dirichlet operator for the lower half-space and has similar representations T - (T) = 
-si 2 m-\I + Sn(T)) = -(1 + s 22(T_))~ 1 s 21 (T).) ' □ 

6. BLOCK TRIANGULAR MATRICES 

So far, everything holds for arbitrary (t independent) coefficients and can be used 

a(x) 
c(x) d(x) 

In the PDE language for L, this means that the conormal derivative is proportional 
to the t derivative: d UA = ad t . If B — A, as in the second section, then this is 
equivalent to B being block lower-triangular as well. 

'a'(x) 
c'(x) d'(x 



in full generality. Assume that A is block lower-triangular: A(x) 



Let us write B(x) 



which we identify with the operator of mul- 



tiplication by B. Then B = UBU also has the same structure 



a 

7 5 



, where a, 7, 5 



are the following operators: if n = 1, they are the operators of multiplication by 
a', c', d! respectively. If n > 2, then a = a', 7 = 1Z1Z*c! and 5 = 1Z1Z*d'lZ7l* . 

Lemma 6.1. If the coefficients are block lower-triangular then S\ 2 (T) is invertible 
from onto H° ± . 

Proof. We write the proof when n > 2. The proof when n = 1 is similar and 
simpler as we do not have to consider the Riesz transforms. By Lemma [5.2[ Si 2 (T_) 

is invertible from T-i^ 1 ^ 2 onto . By the characterization of "H^ 2 using TZ, this is 
equivalent to si 2 (T_)1Z invertible from H -1 / 2 (M Tl ) onto itself. Remark that the same 
holds for S\ 2 {T)1Z from H 1 ' 2 (R n ) onto itself. Thus, looking at the upper off-diagonal 
coefficient in the relation Z3sgn(T)II = sgn(T)B between bounded operators for the 
L? topology, we find 

as l2 {T)mi* = s 12 {T)5. 
Hence we have the equality in £(L 2 (IR n )) 

s 12 (T)K = a'-\s l2 {T)1l){1l*d"R) 
and this implies that si 2 (T_)1Z extends to an invertible operator from (JZ* d'TZ)^ 1 H 1 ^ 2 (J 



onto a'- 1 F 1 / 2 ( 



and 



Using the complex interpolation equalities 

[a'- 1 H l/2 (W n ), H' l/2 (W n )] 1/2 = L 2 (W n ) 

),H- 1/2 (M: n )} 1/2 



\(jvdHiy x H xt2 (} 

we obtain the desired conclusion. 

It remains to explain the interpolation equalities. The first one was actually first 
observed (without proof) in [10] as a consequence of the solution of the 1-dimensional 
Kato square root problem [8] using that a! is an accretive function. See [TTJ Section 
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6] for some account. It can also be seen as a consequence of the Tb theorem. See 
[HI Theoreme 9] for an explicit statement. 

The second equality is a consequence of the solution of the n- dimensional Kato 
square root problem [I] and the results in [5] as follows. Remark that lZ*d'lZ is a 
strictly accretive and bounded operator on L 2 (IR n ). Combining (5j Theorems 5.3, 
7.2, 7.3] using that H s (W l ) is defined by the semi-norms || (— A) S//2 /||2, we have 

[(K*d'Ky l H ll2 (W l ),H- 1,2 {W l )} l , 2 = H° r , 

where T is the sectorial operator (-A) l/2 {lZ*d'TZ). But ii\ = L 2 (R n ) means that 
T has a bounded holomorphic functional calculus on L 2 (R n ), which is the same as 
{lZ*d'lZ){— A) 1//2 has a bounded holomorphic functional calculus on L 2 (IR n ). By [3, 
Theorem 10.1], this is equivalent to the claim 

||l! /2 /I| 2 ^II(-A) 1/2 /I|2, /eDfift 

where Ly is the operator (— A) l / 2 {lZ*d'lZ){— A) 1 / 2 . As Ly is nothing but the diver- 
gence operator — divd'V, the claim is proved in [I]. □ 

Theorem 6.2. The Neumann problem for operators L with block lower-triangular, 
t-independent coefficients A is well-posed for L 2 data. 

Proof. From the relation r^ D = Si^T)" 1 ^ — sn(T)) and the previous lemma, we 
have that T^ D is bounded from 'H L to 'H?.. Thus, given a Neumann data / G L 2 (W l ), 



T^ D f G L 2 (R n , C n ), and from this it is easy to conclude that H { 



o 



/ 



Lu = 



with || .Ho II 2 ~ II /II 2- By Theorem 14.41 (i)> we can define a weak solution u of 
with ||jV(Vu)|| 2 « ||^o II 2 by S7 A u = e-^H . □ 

Remark 6.3. Under the same assumptions, by interpolation, T^ D is bounded from 
U~ s to U~ s for < s < 1/2. When < s < 1/2, given any / G if~ s (M n ) 
the above formalism furnishes a solution to Lw = with conormal derivative / 
and ff R i+nt~ 1+2s \Vt, x u\ 2 dtdx ~ llfll^^s^y By Sneiberg's perturbation result [IB] . 

this can be pushed to s < 1/2 + e for some e > because sn(T) G C{l-i^_ s ) and 
s i 2 (X!) ^ ^(^iT S 5^Z s ) f° r an s e [0) 1]) the spaces are complex interpolation scales 
and invertibility holds at s — 1/2. 

The same arguments apply to block upper-triangular coefficients: this time, the 
lower off-diagonal block coefficient d vanishes but not the upper off-diagonal block 
coefficient b'. 

Theorem 6.4. The regularity problem for operators L with block upper-triangular 
t-independent coefficients A is well-posed for L 2 data. 

Proof. It follows from Theorem 13.41 and Lemma |5~31 that = s 2 i(T) _1 (J — s 22 (T)) 
in £("H^ 2 , l-V] 1 ^ 2 ) from the invertibility of s 2 i(T) from 'H^ 1 ^ 2 onto 'H^ 2 . By 
Theorem 14.41 (i), it is enough to show the invertibility of s 2 i(T) from H° onto 
that is the invertibility of 7?.*s 2 i(T) on L 2 (R n ). This is based on the interpolation 
argument and the formula 

K*s 21 (T)a = (R*d'K)(K*s 2 i(T)) 

in £(L 2 (IR n )) that can be checked as above using that the coefficient d vanishes. 
Details are left to the reader. □ 
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with 7 a M n -valued, measurable and bounded function with div x 7 



Remark 6.5. Observations similar to the ones in Remark 16.31 apply for regularity 
problems with data (the tangential gradient at the boundary) in 'H^ s for block 
upper-triangular coefficients. 

Theorem 6.6. The Dirichlet problem for operators L with block lower-triangular 
t-independent operators A is well-posed for L 2 data. 

Proof. Here, well-posedness is within the class of solutions with J / R i+n t\V t)X u\ 2 dtdx < 

oo and by [T], this implies then that ||A/"(w)|| 2 « JJ R i+n t\ V t)X u\ 2 dtdx ~ ||w|t=o||2- 

Well-posedness follows from the duality principle [2] Proposition 17.6] that the reg- 
ularity problem for L with data in L 2 is well-posed in the class with modified 
non-tangential estimate if and only if the Dirichlet problem for L* with L? data 
is well-posed in the above class. □ 

Remark 6.7. We point out that the relation between L and A is a correspondence 
only when L and A are imposed to be self-adjoint. Otherwise, one L may be rep- 
resented by many different strictly accretive matrices A. As a consequence, there 
are as many Neumann problems for L as choices of A since the conormal derivative 
depends on A. However, for the regularity and Dirichlet problems, the choice of A 
is irrelevant. Although the solution algorithms depend on the coefficients A, they 
produce in the end the same solutions (note that the tangential part of the conormal 
gradient is independent of the coefficients). 

For example, one does not modify L by adding to A any matrix of the form M 7 = 

7(x) r 
-7(z) 

0. First, the divergence free equation implies — div t >X (A + M 7 )Vt )X = — div^AV^ 
in the sense of forms. Second, the accretivity constant of A + M 7 is the lower bound 
of \{A + M 7 + A* + M*), but M 7 + M* = when 7 is W 1 - valued, so this is the same 
as the accretivity constant of A. This argument also shows that A + M 7 remains 
strictly accretive if 7 is C"-valued provided the norm of its imaginary part is not 
too large. 

For systems of size m, this example has to be modified as follows. M n -valued is 
replaced by iS"-valued where S is the space of real and symmetric m x m matrices. 
The divergence free equation is understood coefficientwise Y^j=i djlj = 0- 

Remark 6.8. Recall that any of these well-posedness results is stable under complex 
perturbation in L°° within the class of ^independent coefficients A by [31 Theorem 
2.2] and pQ (the latter reference extends the solution spaces for the results of the 
former to hold). 
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